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1. 


Introduction 


In  order  to  state  clearly  the  purpose  of  this  paper  and  to  motivate  the 
developments  that  follow,  consider  first  the  manner  of  construction  of  an 
infinitesimal  theory  of  deformation.  In  constructing  such  a  theory,  if  as 
usual  it  is  assumed  that  the  displacement  gradient  is  small ,  then  both  the 
strain  and  the  (local)  rotation  are  small  also.  On  the  other  hand,  if  only 
the  strain  tensor  is  assumed  to  be  small,  then  the  rotation  tensor  is  not 
necessarily  small.  It  is  then  of  interest  to  ask  whether  or  not  any  con¬ 
dition  imposed  on  the  strain  field  would  suffice  to  ensure  the  smallness 
of  the  rotation  also.  More  generally,  suppose  it  is  desired  to  have  the 
rotation  moderately  large  in  some  sense.  Is  it  then  possible  to  restrict 
the  strain  field  so  as  to  ensure  that  it  is  accompanied  by  moderate 
rotation? 

Some  insight  into  the  above  questions  is  provided  by  a  well-known 
result  of  the  infinitesimal  theory  which  states  that  the  gradient  of  the 
infinitesimal  rotation  can  be  expressed  in  terms  of  the  gradient  of  the 
strain  field  (see  for  example  Sokolnikoff,  1956,  p.  27).  It  is  not  dif¬ 
ficult  to  see  that  some  corresponding  result  should  exist  even  when  the 
deformation  is  not  infinitesimal.  Indeed,  in  the  context  of  a  finitely 
deformed  body,  the  deformation  gradient  tensor  F  can  be  expressed  as  a 
product  of  the  (local) rotation  tensor  R  and  the  stretch  tensor  U  which  also 
determines  a  measure  of  strain  such  as  the  relative  Lagrangian  strain  E. 

That  there  must  exist  some  connection  between  the  gradient  of  the  stretch 
U  (or  the  gradient  of  the  strain  E)  and  the  gradient  of  the  rotation  R 
becomes  evident  when  one  recalls  the  compatibility  condition  that  F  must 
satisfy.  The  main  purpose  of  the  present  paper  is  to  derive  a  representa¬ 
tion  for  the  rotation  and  deformation  field  directly  from  the  strain,  which 
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can  then  be  used  to  estimate  the  magnitude  of  R  in  terms  of  the  magnitude 
of  U  or  E.  Results  of  this  kind  are  •'f  interest  in  various  contexts  and 
especially  for  materials  undergoing  deformation  in  which  the  strain  and 
rotation  are  not  necessarily  of  the  same  orders  of  magnitude.  In  fact,  in 
a  number  of  theories  for  special  bodies,  such  as  those  for  shells  and 
rods,  there  are  circumstances  in  which  the  motion  results  in  small  strain, 
while  the  rotation  and  the  deformation  may  be  large  or  moderately  large. 

Motivated  in  part  by  the  remarks  made  in  the  preceding  paragraph, 
most  of  the  paper  is  concerned  with  the  development  of  a  procedure  for  the 
construction  of  a  theory  of  infinitesimal  strain  accompanied  by  moderate 
rotation,  although  the  procedure  is  applicable  to  other  situations  in  which 
the  strain  and  the  rotation  may  be  of  different  orders.  Thus,  after 
introducing  a  geometrically  appealing  definition  for  moderate  rotation 
and  establishing  a  number  of  relevant  theorems  pertaining  to.  its  properties, 
estimates  are  obtained  for  the  angle  of  rotation  and  the  moderate  rotation 
tensor  and  these  are  eventually  expressed  in  terms  of  infinitesimal  strain. 
These  and  related  kinematical  results  are  then  put  in  a  properly  invariant 
form,  and  the  invariance  of  constitutive  equations  under  arbitrary  superposed 
rigid  body  motions  is  discussed  in  the  case  of  an  elastic  material  under¬ 
going  small  strain  accompanied  by  moderate  rotation.  Although  throughout 
the  paper  use  is  made  of  the  direct  (coordinate- free)  notation,  which  often 
allows  the  results  to  be  stated  in  their  simplest  form,  on  occasions  we 
also  employ  the  component  forms  (at  least  partially)  of  the  various  equa¬ 
tions  since  these  are  more  convenient  in  explicit  calculations.  A  brief 
account  of  the  notations  used  and  some  mathematical  preliminaries  is  given 
at  the  end  of  this  section  and  additional  notations  and  mathematical 
terminology  are  collected  in  Appendix  A. 
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The  problem  of  the  determination  of  the  rotation  tensor  R  from  c 
strain  tensor,  and  the  connection  between  them,  has  a  long  history  that 
dates  back  o  Cauchy.  A  detailed  account  of  the  subject  can  be  found  in 
sections  34-38  and  55-57  of  7i*uesdell  and  Toupin  (1960).  Truesdell  and 
Toupin  (1960,  p,  276)  also  discuss  in  some  detail  a  measure  of  rotation  -- 
called  mean  rotation  --  introduced  by  Novozhilov  (1953)  and  use  this  in 
their  development  of  a  theory  of  infinitesimal  strain  and  infinitesimal 
rotation  (Truesdell  and  Toupin  1960,  p.  305).  More  recently,  in  the  con¬ 
text  of  finite  deformation,  the  problem  of  the  determination  of  the  rotation 
tensor  R  from  the  strain  tensor  has  been  discussed  by  John  (1961)  and  by 
Shield  (1973).  In  particular,  Shield  has  derived  a  useful  representation 
of  the  integrability  relations  for  finite  strain  in  a  Euclidean  space,  and 
has  further  shown  that  for  two-dimensional  deformation  the  rotation  and 
deformation  field  can  be  determined  directly  from  the  knowledge  of  strain 
by  a  line  integral. 

1 . 1  Scope  and  outline  of  contents 

After  collecting  some  kinematical  and  kinetical  results  in  section  2, 
for  clarity  and  ease  of  reference  we  include  a  sketch  of  the  derivation  of 
the  differential  equations  for  rotation  R  in  section  3.  Since  the  presenta¬ 
tion  of  this  development  differs  somewhat  from  that  given  by  Shield  (1973), 
one  or  two  of  the  details  are  collected  in  the  first  part  of  Appendix  B 
(between  (Bl)  and  (B4))  whore  a  correspondence  with  Shield's  (1973)  main 
results  are  indicated.  We  also  take  this  opportunity  to  include  in 
section  3  two  elementary  kinematical  theorems  concerning  the  rotation  R 
and  the  deformation  gradient  which  are  utilized  later  in  the  paper. 

We  begin  section  4  by  defining  in  precise  terms  four  measures  of 
smallness  denoted  by  and  associated,  respectively,  with 


strain,  strain  gradient,  rotation  and  rotation  gradient  [see  the  expressions 
(4.1) -(4. 2)  and  (4.5) - (4 .6) ] .  We  then  observe  the  conditions  under  which  these 
four  measures  can  be  related  [see  (4.7)]  and  go  on  to  show  that  the  rotation 
gradient  is  of  the  same  order  of  magnitude  as  the  strain  gradient,  i.e., 
is  of  the  same  order  as  [see  (4.17)].  Under  the  additional  assumption 
that  the  rotation  tensor  is  the  identity  tensor  at  some  material  point  in 
the  body,  we  prove  in  subsection  4.1  that  the  rotation  tensor  is  of  the 
uame  order  of  magnitude  as  the  strain  gradient,  i.e.,  e2  is  of  the  same 
order  as  [see  (4.26)].  Later,  making  use  of  a  further  assumption  that 
relates  the  two  measures  of  smallness  eQ  and  £^,  we  express  the  rotation 
tensor  and  all  other  kinematical  quantities  in  terms  of  or  equivalently 
the  infinitesimal  strain. 

In  the  second  part  of  section  4,  we  introduce  two  separate  definitions: 

One  for  (1)  infinitesimal  strain  accompanied  by  infinitesimal  rotation 
(Definition  4.1)  and  another  (2)  for  infinitesimal  strain  accompanied  by 
moderate  rotation  (see  Definition  4.2).  Although  our  main  objective  here 
is  to  deal  with  the  latter,  we  find  it  instructive  to  prove  a  number  of 
results  in  subsection  4.2  pertaining  to  the  infinitesimal  theory  of  motion. 

The  remainder  of  section  4  (subsections  4.3  and  4.4)  deals  entirely  with 
the  case  of  infinitesimal  strain  accompanied  by  moderate  rotation.  In 
Theorem  4.3,  we  essentially  prove  that  if  the  strain  is  of  0(eQ)  and  R  is 

Tx 

a  moderate  rotation,  then  the  relative  displacement  gradient  is  of  0(ep. 

In  a  subsequent  Theorem  4.4,  with  the  help  of  a  known  representation  for  R 
in  terms  of  angle  of  rotation  9  and  the  results  stated  in  the  preceding 
paragraph,  we  estimate  the  order  of  magnitude  of  R  in  terms  of  the 
infinitesimal  strain  We  also  obtain  in  subsection  4.3  expressions  for 
the  deformation  gradient  F  and  the  relative  displacement  gradient  H  when 
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the  motion  is  such  that  the  strain  is  of  Q(eo)  but  is  accompanied  by 

% 

moderate  rotation  of  0(eQ) . 

In  a  recent  paper,  Casey  and  Naghdi  (1981)  have  constructed  a 
properly  invariant  theory  in  which,  apart  from  superposed  rigid  body 
motions,  both  the  strain  and  rotation  (and  hence  also  the  relative 
displacement  gradient)  are  infinitesimal  and  have  further  shown  that  in 
this  infinitesimal  theory  of  motion  the  constitutive  equations,  as  well 
as  the  equations  of  motion,  are  all  properly  invariant  under  arbitrary 
(not  necessarily  infinitesimal)  superposed  rigid  body  motions.  In  the 
latter  part  of  section  4  (see  subsection  4.4),  by  using  the  procedure  of 
Casey  and  Naghdi  (1981)  we  remove  from  all  motions  the  translation  and 
rotation  at  any  particle  Y  of  the  body,  thereby  rendering  all  kinematical 
quantities  (including  the  infinitesimal  strain)  properly  invariant  under 
superposed  rigid  body  motions.  Finally,  in  section  5,  in  a  manner  similar 
to  that  of  Casey  and  Naghdi  (1981)  we  briefly  discuss  the  invariance 
properties  of  the  constitutive  equations  for  small  strain  accompanied  by 
moderate  rotation  and  further  demonstrate  in  Appendix  C  that,  as  in  the 
case  of  the  infinitesimal  theory,  for  the  theory  under  discussion  the 
particular  choice  of  the  particle  Y  is  again  immaterial  in  effecting  the 
properly  invariant  nature  of  the  results. 

1 . 2  Notation  and  mathematical  preliminaries 

We  close  this  section  with  a  short  glossary  of  notations  and  some 
mathematical  terminology  that  will  be  needed  in  what  follows.  Any  linear 
mapping  from  V,  the  three-dimensional  translation  vector  space  associated 
with  the  Euclidean  point  space  £,  into  V  will  be  called  a  second  order 
tensor.  The  trace  and  determinant  functions  of  second  order  tensors  will 
be  denoted,  respectively,  b>  tr  and  det.  The  transpose  of  a  second  order 
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T 

tensor  T  will  be  denoted  by  T  ,  while  the  inverse  of  T  if  it  exists  will 
be  denoted  by  T~*.  The  usual  inner  product  on  V  is  written  a*  b  for  any 
two  vectors  a,b  £  V  and  the  (induced)  norm,  or  magnitude,  of  a  is  given  by 

Is 

||a[j  =  (a*  a)*.  The  set  of  second  order  tensors  can  be  provided  with  an 
inner  product  A*  B  =  tr(Al8)  and  a  norm  ||a||  *  (A  •  A) 2  for  any  second  order 
tensors  A  and  B.  We  note  that  the  definition  of  the  norm  of  a  second  order 
tensor  satisfies  the  usual  properties  of  the  norm,  i.e., 

II  a|5  >  0  ,  and  ||  A||  =0  if  and  only  if  A  =  0  , 

IMi  =  |ct|  || a||  ,  ci.i) 

l|A  +  B||  <  11  All  *  U  B I) 

for  any  second  order  tensors  A  and  B  and  for  any  scalar  a.  Further,  it 

*«»  «** 

can  be  shown  that 

l|A  ai|  <  I)  A |l  U a ||  (1.2) 

for  any  second  order  tensor  A  and  any  vector  a.  In  addition,  for  any 
second  order  tensors  A  and  B,  we  have 

II A  B||  <  ||A||||B||  .  (1.3) 

Ihe  tensor  product  a®b  of  any  two  vectors  a,b  €  v  is  the  second  order 
tensor  defined  by  (a®b)v=  (b  •  v)a  for  every  vector  v.  We  recall  the 
formulas  tr (a ® b)  =  a  •  b,  (a®b)"^  =  b®a,  (a®b)(c®d)  =  (a®d)(b,c), 

(a®b)  •  (c®d)  =  a  •  c  b  •  d,  which  hold  for  all  vectors  a,b,c,d  €  V.  In 
addition,  we  note  that 

||a®b||  =  Hall  SlbH  (1.4) 

for  all  a,b  €  v. 
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2.  Some  preliminary  kinematical  and  kinetical  results 

Let  X  be  a  particle  of  a  body  iB  and  denote  by  X  and  x  the  position 
vectors  of  X  in  a  fixed  reference  configuration  kq  and  in  the  current  con¬ 
figuration  k,  respectively.  A  motion  of  SB  is  the  mapping  x  which  assigns 

the  position  vector  x  =  x(X»t)  to  the  particle  X  at  time  t.  In  the  con- 

* 

figuration  kq,  let  SB  occupy  a  region  SRo  embedded  in  a  3-dimensional 
Euclidean  space  £  and  denote  the  boundary  of  SRq  by  3<Rq.  The  image  of 
SRq  in  <  will  be  denoted  by  SR  with  boundary  3 SR.  We  asstime  that  at  each 
fixed  t,  the  mapping  x  of  into  SR  possesses  a  smooth  inverse  denoted 
by  X*1. 

*v  < 

X 

Let' F= Grad  x  be  the  deformation  gradient  relative  to  the  configura¬ 
tion  Kq  and  recall  that  det  F>0.  Then,  by  the  polar  decomposition  theorem, 

F  can  be  decomposed  in  the  form 

F  =  R  U  ,  (2.1) 

where  the  right  stretch  U  is  a  symmetric  positive  definite  second  order 
tensor,  the  (local)  rotation  R  is  a  proper  orthogonal  tensor  satisfying 

r‘r  =  R  RT  =  I  ,  det  R  =  1  (2.2) 

and  I  stands  for  the  identity  tensor.  The  right  Cauchy-Green  measure  of 
deformation  C  can  be  represented  in  the  forms 

C  =  FTF  =  U2  (2.3) 

and  the  relative  (Lagrangian)  finite  strain  tensor  E  is  given  by 

E  =  %(C-I)  .  (2.4) 

* 

A  region  is  regarded  here  as  a  nonempty  connected  and  compact  subset  of  £ 
having  a  piecewise  smooth  boundary. 

Definition  of  the  gradient  of  a  vector  function  is  given  by  (A4)^  of 
Appendix  A. 
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In  order  to  express  certain  expressions  in  component  forr.,  it  is 


for  some  proper  orthogonal  second  order  tensor-valued  function  Q(t)  of  time, 
some  vector-valued  function  a(tl  of  time  and  some  real  constant  a.  The 
configuration  of  S,  at  time  t+  in  the  motion  x+  will  be  denoted  by  k+. 

The  deformation  gradient  F+(X,t+)  calculated  from  (2.9)  is  related  to  F  by 

F+  =  Q(t)F  . 

Then,  it  can  be  easily  verified  from  (2.3),  (2.4)  and  (2.1)  that 

=  c  ,  U+  =  U  ,  E+  =  E  ,  E^k  =  E  K  ,  R*  =  Q(t)R  . 

The  relative  displacement  field  associated  with  x+  is  u+  =  x+  -  X  and  its 

■4*  4* 

gradient  H  = Grad  u  in  terms  of  (2.6)  is  given  by 

H+  =  QF  -  I  =  QH  +  Q  -  I  . 

Thus,  C,U,E  are  unaltered  under  superposed  rigid  body  motions,  while  R  is 

§ 

unaltered  apart  from  orientation  .  The  displacement  gradient  H,  however, 
is  neither  unaltered  nor  unaltered  apart  from  orientation  under  the 
transformation  (2.9). 

For  later  reference,  we  recall  here  an  exact  representation  for  the 

★ 

rotation  tensor  R,  namely 

R  =  I  +  4>  +  T  , 

where 

^For  the  motivation  and  the  precise  meaning  of  this  terminology,  see 
Green  and  Naghdi  (1979) . 

•b 

oee  Truesdell  and  Toupin  (1960,  p.  20,  Eq.  (37.17)). 
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(2.10) 


(2.11) 


(2.12) 


(2.13) 


2  =  2*  =  (1  -  COS  9)  (y  ®  y  -  I)  =  (1  -  cos  9)  (yRyj  -  6KL)eK  ®  <?L  , 

!  =-f  ="sin  G  £  H  =-sin  9  £klmum!k  ®  ?l  • 

In  (2.14)^  2«  <5^  are  t^e  coraPonents  of  the  identity  tensor  I, 

e  =  e^MeK  ®  ®  eM  is  tiie  Permut.ation  symbol  in  3-space,  y  =  y^eM  is 

the  unit  vector  in  the  principal  direction  of  R  that  is  associated  with 
the  principal  value  1,  and  9  lepresents  the  angle  of  rotation  calculated 
from  the  relationship 

tr  R  =  1  +  2  cos  9  ,  -tt  <  9  <rr  . 


We  introduce  now  some  kinetical  quantities  needed  in  the  development 
of  section  5-  Thus,  with  reference  to  the  motion  x>  let  p  =  p(X,t)  be  the 
mass  density  in  the  configuration  tc,  n  the  outward  unit  normal  to  the 
surface  3<R,  t  =  t(X,t  ;  n)  the  stress  vector  acting  on  this  surface  and 
T  =  T(X,t)  the  associated  Cauchy  stress  tensor.  The  corresponding  quantities 
in  the  configuration  k+  will  be  denoted  by  p+,n+,t+  and  T+,  respectively. 

We  recall  that  under  the  superposed  rigid  body  motions  (2.9) ,  p+  and  n+ 
transform  according  to 


(2.14) 


(2.15) 


p+  =  p  ,  n+  =  Q(t)n 


and  adopt  the  usual  assumption  that  t+  for  he  motion  x  is  related  to  t 


t  =  Q(t)t  . 


It  then  follows  from  the  relation  t  =  T  n,  (2.16)^  and  (2.17)  that 

T+  =  Q(t)T  QT(t)  . 


(2.16) 


(2.17) 


(2.18) 


■■  .  ‘~tneaffg.it 


3. 


Differential  equations  for  the  rotation  tensor 

In  order  to  examine  locally  the  deformation  of  the  body  S,  consider  a 

material  line  element  dX  at  X  in  the  reference  configuration  <  .  As  a 

~  ~  ~o 

result  of  deformation,  dX  is  mapped  into  a  line  element  dx =  F  dX  at  x  in 
the  present  configuration  k.  Since  the  space  C  is  3-dimensional,  in  order 
to  determine  the  deformation  of  an  arbitrary  material  line  element  it  will 
suffice  to  consider  the  deformation  of  three  mutually  orthogonal  line 
elements  dXM,  namely 


d^  =  F  d^  ,  (M-  1,2,3)  .  (3.1) 

Let  the  magnitudes  of  dX^  and  dx^  be  denoted  by  dS^  and 

ds^,  respectively,  and  introduce  the  unit  vectors  in  the  directions  of 

dX^  and  the  unit  vectors  k^  in  the  direction  of  dx„.  Hence,  the  three 

unit  vectors  (M  =  1,2,3)  form  an  orthonormal  basis  which,  without  loss 

in  generality,  will  be  identified  with  the  orthonormal  basis  {e,,}  in  <  . 

~o 

In  general,  the  line  elements  dX^  undergo  both  stretch  and  rotation  and 
the  ratios  (ds^/dS^  ds2/dS2,  ds3/dS3)  denoted  by  (M =  1,2,3)  are  called 
the  stretch  of  the  line  elements.  The  above  observations  may  be  summarized 
as 

dsM 

d?M  =  ?M  dSM  *  dhi  =  *Lm  dsM  ’  XM  =  dS“  Cn°  Sum  °n  M)  ■  (3>2) 

M 

It  follows  from  (3.1)  and  (3,2) ^  2  ^  that 

d?M‘  !i  =  FiM  dSM  =  XM  kMi  dSM  (n0  sum  on  M)  *  (3,3) 

where  k^  are  the  components  of  the  unit  vector  k^  referred  to  the 
orthonormal  basis  {e^}  in  the  current  configuration  <  and  F^  are  the 
components  of  F  referred  to  the  basis  {e^  ®  e^} ,  i.e., 


F.„e .  ® eu  =  e. 
iM-i  ~M  -M 


f-Sl  =  FiM~i 


(3.4) 
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In  view  of  (3.4),  the  results  (3.3)  nay  also  be  displayed  as 


$ki  =  (no  sum  on  M)  . 


(3.S) 


It  should  be  noted  that  since  F  is  invertible,  are  linearly  independent 

and  form  a  basis.  Moreover,  Mj^  can  be  interpreted  as  the  image  of  the 

material  basis  (or  e^)  in  the  configuration  k.  It  is  clear  that  the  knowledge 

of  the  images  of  IC|  at  every  point  gives  full  information  about  the 

components  ^vi =  ei  *  ^  e"  55  (or  t^e  comPonents  F^M) ,  the  necessary 

and  sufficient  conditions  for  the  existence  of  a  deformation  function 

* 

whose  gradient  is  F  are  the  compatibility  conditions 

T1 

(Grad  F)  =  Grad  F  or  F.„  r  =>  F..  „  .  (3.6) 

—  ~  1  ,  Lj  1  L  ,  N 

In  terms  of  M. the  conditions  (3.6)  can  also  be  written  as 


(3.7) 


Clearly  the  inner  product  of  and  gives  the  Cauchy-Green  measu-e,  i.e., 


C„.  =  F.  ,,F. .  =  F.„e.  •  F . .  e  .  = 

KL  lK  lL  iK~i  ~K 


M, 


(3.8) 


and  it  is  seen  that  the  components  of  the  metric  tensor  associated  with  the 
material  basis  are  C^. 

Given  C^,  it  does  not  necessarily  follow  that  there  exists  a  deforma¬ 
tion  function  x  whose  gradient  F  satisfies  (2.3)^.  Thus,  we  must  ask  the 
following  question:  If  C  is  specified  at  every  point  of  the  body,  can  we 
find  F  (or  equivalently  M^)  derivable  from  x  such  that  (2.3)1  or  (3.8) 
are  satisfied?  To  this  end,  we  proceed  to  express  the  derivative  of  in 
terms  of  the  derivatives  of  C.  First,  we  differentiate  (3.8)  with  respect 
to  and  obtain 

The  definition  of  the  transpose  in  (3.6) ^  is  given  by  (A2)^  of  Appendix  A. 
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CKL,M  "  (5?K  *  ~L5  ,M  =  ~K,M  *  ~L  +  ~K  *  ~L, 


M  ' 


(3.9) 


Next,  by  suitable  combination  of  expressions  of  the  type  (3.9)  and  with  the 
help  of  (3.7),  the  derivative  of  Mj,  can  be  expressed  in  terms  of  partial 
derivatives  of  components  alone,  i.e., 

~L*~K,M  “  FiLFiK,M  =  ^CCLM,K  +  CLK,M' CKM,J  (3,10) 

d*J* 

or  equivalently  in  coordinate- free  notation1' 

T  T 

FT  Grad  F  *  [(Grad  C)  1  +  Grad  C-  (Grad  C)  3]  .  (3.11) 


Thus,  if  there  exists  an  F  which  satisfies  (3.6)  and  (2.3)^,  then  F  must 
satisfy  the  differential  equations  (3.10)  or  (3.11).  Clearly,  the  dif¬ 
ferential  equations  (3.10)  or  (3.11)  are  the  necessary  conditions  for  the 

determination  of  F  from  the  knowledge  of  C.  The  converse,  proved  in 
* 

Appendix  B,  is  :  if  a  field  F  satisfies  (3.10),  then  (2.3)^  holds  and  further¬ 
more  a  motion  x  exists  whose  gradient  is  F. 

Since  the  stretch  tensor  U  is  uniquely  determined  by  C,  we  may  use 

s 

(2.1)  to  obtain  the  following  differential  equations  for  the  rotation  R 
in  terms  of  U: 


R1  Grad  R  =  J 


(3.12) 


where  J  is  a  third  order  tensor  given  by 


.1 


§ 


The  notations  Ti,T2  and  Tj  for  the  transpose  of  third  order  tensors  is 
introduced  in  (A2)  of  Appendix  A. 

k 

This  sufficiency  argument  is  discussed  in  the  first  paragraph  of  Appendix  B. 

A  sketch  of  the  derivation  of  the  differential  equations  (3.12)  or  equiva¬ 
lently  (3.14)  is  given  in  the  second  paragraph  of  Appendix  B  leading  to  (B4) 
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i— itV"  mcza— a . Trrra  w 


j-.  j  a*.  t,  . .  - 


J[  =  -Jl  2  =  ^{if1  [Grad  U-  (Grad  U)  3] 

T  T 

+  ([Grad  U-  (Grad  U)  ^U'1)  1 

•v  ■»*  — 

T  T 

+  (U_1  [Grad  U  -  (Grad  U)  3]U_1)  2U}  .  (3.13) 

•V  **  «v  *** 

Thus,  if  R  exists  such  that  (2.1)  holds,  then  it  is  necessary  that  R  satisfy 

(3.12) .  Conversely,  if  a  proper  orthogonal  tensor  R  is  a  solution  of  (3.12), 
then  by  reversing  the  procedure  used  in  obtaining  (3.12)  from  (3.10),  it  can 
be  shown  that  the  tensor  F=R  U  satisfies  (3.10)  and  the  existence  of  a 
deformation  function  x  is  ensured. 

In  the  remainder  of  the  paper,  it  is  often  more  convenient  to  write 

(3.12)  as 

RTR.k  ■  Ak  ,  (K- 1,2,3)  ,  (3.14) 

where  each  (K  *  1,2.3)  is  a  skew- symmetric  second  order  tensor  such  that 

~K  =  "  6k  =  ^  !k 

=  !i{U_1U(K  -  U  KU_1  +  U(L(eK  ®  eL)U_1 

-  U’ 1  (e.  9  eJU  .  +  U_1U  rU(e„  ©  e.)lf 1 

-  U~ 1  (e.  ©  eJU  U  .if1}  .  (3. IS) 

•v  ■** 1\  *»  ‘“jL'- 

It  should  be  noted  that  the  right-hand  side  of  (3.13),  or  equivalently 
(3.15),  involves  only  the  stretch  U  and  its  gradient. 

Integrability  conditions  for  the  existence  of  a  solution  R  to  (3.14) 
have  been  discussed  by  Shield  (1973)  based  on  a  paper  of  Thomas  (1934).  A 
brief  statement  of  these  conditions  in  the  context  of  the  present  paper  is 
given  following  (B4)  of  Appendix  B. 

In  the  rest  of  this  section,  we  state  two  theorems,  the  proofs  of 
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which  utilize  the  representation  (3.14)  or  some  of  the  preceding  results. 

The  first  of  these  pertains  to  uniqueness  of  the  rotation  tensor  to 
within  a  rigid  body  rotation  and  may  be  stated  as 

TheoremSjl.  Let  a  stretch  field  U  derivable  from  a  motion  X  exist 
such  that  the  compatibility  conditions  (3.6)  are  satisfied  and  let  and 
R~  be  any  two  rotation  tensors  satisfying  (2,2).  Then,  R,  and  R«  cor- 
responding  to  U  are  solutions  of  (3.14)  if  and  only  if  they  differ  by  a 
rigid  body  rotation. 

Proof.  We  first  prove  the  necessity.  Suppose  that  R.  and  R_  are 
any  two  distinct  proper  orthogonal  tensors,  each  of  which  satisfy  (3.14) 
corresponding  to  the  same  stretch,  i.e., 

?1?1,K  *  6k  *  ?2?2,K  “  •  (3,16) 

T 

Consider  next  the  space  derivative  of  the  product  R-R-,  namely 

- 1 

<5i5J>,k  ’  5Jc5l51,K>sJ*Si<5l.KS2)£l 

-  6k~2  '  ?1^K?2 

-  0  ,  (3.17) 

where  in  obtaining  (3.17)^  the  first  term  has  been  premultiplied  by 

T  T 

RjRj  *  I  and  the  second  term  has  been  postmultiplied  by  ^2R2 *  *  where 

use  has  been  made  of  (3.16)^  £  and  (3.15)^.  It  follows  from  (3.17)  that 
T 

the  product  must  be  a  second  order  tensor  which  is  independent  of 

position  but  may  be  a  function  of  time.  Since  each  of  the  two  tensors  R^ 

T 

and  R-,  is  proper  orthogonal,  the  product  R^R^  is  also  proper  orthogonal, 

T 

say  Q0(t).  Thus,  we  may  conclude  that  R^R,  «Qq  or  equivalently 

h  "  9o*2  *  (3-18) 

15. 


Hence  R,  and  R_  differ  at  most  by  a  proper  orthogonal  tensor  function  of 
time  corresponding  to  a  superposed  rigid  body  rotation.  To  show  that  the 
conclusion  reached  is  also  sufficient,  we  only  need  to  replace  R  in  (3.14) 
by  Q  R  and  the  theorem  is  proved. 

An  immediate  consequence  of  Theorem  3.1  is  the  following 

CorollaryJLl .  Corresponding  to  a  given  field  of  deformation  C,  the 
deformation  gradient  field  F  can  be  determined  to  within  a  rigid  body 
rotation. 

The  proof  follows  at  once  by  recalling  (2.1)  and  noting  that  corresponding 
to  a  given  C  the  tensor  U  is  uniquely  determined  by  (2.3), •  A  proof  of 
this  result  was  apparently  first  given  by  Shield  (1973,  p.  484)  who  employs 
a  different  procedure  than  that  of  our  Theorem  3.1  and  does  not  make  use 
of  (3.14). 

In  the  statement  of  Theorem  3.1,  the  existence  of  a  motion  x  whoso 
stretch  field  is  U  was  assumed.  Suppose  instead  that  a  field  U  is  prescribed 
as  a  function  of  X.  It  is  then  natural  to  ask  under  what  conditions  such 
a  stretch  is  the  gradient  of  some  deformation  function  x*  In  other  words, 
what  restrictions  must  be  placed  on  1)  in  order  to  guarantee  that  the 
deformation  gradient  is  of  the  form 


r  ***- 


Proof,  Suppose  a  motion  exists  whose  gradient  satisfies  (3.19). 

Then,  the  compatibility  conditions  (3.6)  imply  the  restrictions  (3.20)  as 

necessary  conditions  for  the  existence  of  x*  To  show  sufficiency,  assume 

the  restrictions  (3.20)  on  U.  With  the  use  of  (3.20),  it  can  then  be  shown 

★ 

that  the  right-hand  side  of  (3.15)  vanishes  identically  and  hence  A^-0. 

T 

Further,  by  (3.14)  we  have  R  R  »  0.  This,  in  turn,  implies  that  R  ,.*0 
since  R  is  nonsingular.  Hence,  R  must  be  a  function  of  time  only  cor¬ 
responding  to  a  superposed  rigid  body  rotation,  say  Qq.  Thus,  a  motion 
T  T 

Q  y  whose  gradient  is  OF  is  sufficient  to  satisfy  (3.19)  and  the 
theorem  is  proved. 
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Kinematics  of  deformation  with  moderate  local  rotation 


We  are  concerned  here  with  a  kinematical  development  in  which  the 
Lagrangian  strain  is  infinitesimal  but  is  not  necessarily  accompanied  by 
infinitesimal  rotation.  To  make  these  notions  precise,  in  a  manner 

similar  to  that  of  Casey  and  Naghdi  (1981),  we  first  define  a  measure  of 

6 

smallness  of  strain  by  the  nonnegative  real  function 

e  =  e  (t)  -  sup  ||E(X,t)||  ,  (4.1) 

0  0  X  €  jR  ~  ~ 

**  0 

where  sup  stands  for  the  supremum  (or  least  upper  bound)  of  a  nonempty 

bounded  set  of  real  numbers.  If  h  (E)  is  any  scalar-,  vector-,  or  tensor- 

valued  function  of  E  defined  in  the  neighborhood  of  E  *  0  and  satisfying 

the  condition  that  there  exists  a  nonnegative  real  constant  C  such  that 

II h  (E)||  <  Ce^  as  e  -*•(),  then  we  write  h  «0(en)  as  e  +0. 

The  statement  E  *  0(e  )  as  £  -*-0  does  not  imply  any  restriction  on  the 
~  -  o  o 

space  or  time  derivatives  of  E.  In  particular,  it  is  possible  for 
Grad  E  =  E  „  ® e,.  or  equivalently  E  v  to  be  finite  while  E  itself  remains 
small.  To  deal  with  such  circumstance: ,  *  introduce  a  second  measure  of 

'It 

smallness  by  a  nonnegative  real  function 

£  1  =  e,  (t)  =  max  sup  jj E  (X, t)  ||  ,  (K  *  1,2,3)  .  (4.2) 

K  x  e  <R  ■*’ 

-  0 

If  h1  (E  v)  is  any  scalar-,  vector-,  or  tensor-valued  function  of  Grad  c 

g 

sIn  the  paper  of  Casey  and  Naghdi  (1981),  a  quantity  corresponding  to 
was  defined  in  terms  of  the  displacement  gradient  rather  than  the  strain. 

★ 

In  writing  (4.2),  for  simplicity  we  have  used  the  norm  of  the  second  order 
tensors  E  ^  defined  in  section  1.  It  is,  of  course,  possible  to  define 
£l  in  terms  of  Grad  E,  but  this  requires  also  a  definition  for  the  norm 
of  a  third  order  tensor  which  is  not  introduced  in  section  1  or 
Appendix  A. 
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defined  in  the  neighborhood  of  E  v  ■  (E  .  ,E  ,,E  _)  *  0  satisfying  the  con- 
dition  that  there  exists  a  positive  real  constant  D  such  that 
|1  h1  (E  ^)  ||  <  De^  as  t^  +  0,  then  we  write  h^Ofe")  as  ^->-0.  Further,  if 
there  exists  a  nonnegative  constant  C"  such  that 

£1  <  ^eo  wit^  *nteSers  *  (4.3) 

then  it  can  be  shown  that 

Hhj. Cg.K3  il  <  tCCn/k)e||A/k  •  (4.4) 

Provided  that  (4.3)  holds,  the  last  result  implies  that  any  function 
h^E  K)  of  O(e^)  as  e^O  is  also  of  0(0^^)  as  eQ  +  0. 

We  now  introduce  two  additional  measures  of  smallness,  one  associated 
with  che  rotation  tensor  and  another  with  its  gradients.  Thus  we  define 
the  measures  of  smallness  e2  and  e^,  respectively,  by  the  nonnegative 
real  functions 

e9  3  £9(t)  3  sup  1|  R  -  1 1|  (4.5) 

Z  Z  X  6  <R0  ~  ~ 

and 

e,  3  e,(t)  =  max{  sup  || R  J|  }  ,  (K- 1,2,3)  .  (4.6) 

o  o  K  X  €  <R  ■',K 

~  o 

If  h_(R-I)  is  any  scalar-,  vector-,  or  tensor-valued  function  of  R  defined 

it  **  “ 

in  the  neighborhood  of  R*  I  satisfying  the  condition  that  there  exists  a 
positive  real  constant  D 2  such  that  ||h2(R-I)||  <  D2e2  as  e2  +  0,  ^ien  we 
write  h2*0(e2)  as  e2  +  0.  Similarly,  if  hj(R  K)  is  any  scalar-,  vector-, 
or  tensor-valued  function  of  R  ^  (K  ■  1,2,3)  defined  in  the  neighborhood 

The  details  are  discussed  following  (BS)  in  Appendix  B. 
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of  R  K  =  0  satisfying  the  condition  that  there  exists  a  positive  real 
constant  D3  such  that  ||h3(R  K)||  <  as  e3  +  0,  then  wc  write 
h3  =  0(e3)  as  e3  +  0.  Further,  the  four  measures  em,  (m=o,l,2,3)  can 
be  related  if  there  exist  nonnegati^ constants  C  ,  m,n  =  o,l,2,3, 
such  that 

£m  <  ^mn£n  with  k,£  inte8ers  *  (4.7) 

and  then  similar  to  (4.4)  it  can  be  shown  that  any  function  of  0(en)  as 

m 

e  •+•  0  is  also  of  OCe11^)  as  e  -*-0. 
m  n  J  n 

Since  R  is  a  proper  orthogonal  tensor,  it  follows  from  the  definition 
of  the  form  of  a  second  order  tensor  that 

II R-I  il 2  “  tr{(R-I)T(R-I)}  =  2(3  -  tr  R)  .  (4.8) 

Recalling  the  expression  for  tr  R  in  (2.15),  which  implies  that  tr  R 

N* 

satisfies  the  inequality  -1  <  tr  R  <  3,  it  follows  from  (4.8)  that 
2 

0  <  ||  R-I  ||  <8  and  consequently  the  measure  of  smallness  e2  in  (4.5) 

is  bounded  from  above  by  2/2",  i.e.,  e2  <  2/2". 

A  solution  of  (3.14)  involves  both  U  and  its  derivatives  U  r  In 
order  to  accommodate  the  strains  and  strain  gradients  of  different  orders 
of  magnitude,  we  first  express  U  ..  in  terms  of  E  and  E  ...  Remebering  that 
by  (2.3)2  and  (2.4)  the  stretch  U  may  be  regarded  as  a  function  of  E,  from 

'it  k 

the  Taylor  expansion  of  U(E)  about  E  =  0  follows 

U  *  I  +  E-  jy  E2+yj-  (3E3)  +  •••  (4.9) 


A  statement  of  Taylor's  formula  of  the  form  (4.9)  is  given  by  Theorem  8.14.3 
of  Dieudonne  (1969,  p.  190).  Since  all  the  spaces  considered  here  are  real 
Euclidean  spaces  and  since  all  finite  dimensional  Euclidean  spaces  are 
Banach  spaces,  the  theorem  in  (DieudonnS  1969)  is  directly  applicable. 
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and  by  differentiation  we  have 


U  v  «  E  „  -  IfiCE  VE  +  EE  „) 

j  N  ~  j  N  *«mj  I  n 


+  Jj[E  „E2  +  EE  „E  +  E2E  J  + 

■W  I  l\*v  p  ^  |  rt 


(4.10) 


where  in  writing  (4.9)  we  have  also  used  the  fact  that  U*  I  at  E  =  0.  Before 
proceeding  further,  in  the  context  of  the  classical  infinitesimal 
kinematics,  we  recall  the  approximate  formulae  for  U,C  and  their  inverses, 
all  estimated  in  terms  of  the  infinitesimal  Lagrangian  strain,  i.e., 


E  ■  0(e  )  ,  U  =  I  +  E  =  l  +  0(e  )  ,  U_1  *  I  -  E  =  I  -  0(e  )  , 

-  ~  O  *s»  —  «*.  0  -M  —**•*+,  0 

C  =  I  +  2E  =  I  +  0(e  )  ,  C_1  =  I  -  2E  =  I  -  0(e  as  e  ■*•0  . 

M  M  >V  0  ~  ~  *»  «*  O  O 

Now  suppose  that  E  =  0(eQ)  as  £Q>0  in  (4.10)  but  as  yet  impose  no  restric¬ 
tion  on  E  It  follows  that  each  term  of  the  last  bracket  in  (4.10)  is 
2  2 

of  0(e  )  and,  to  the  order  e.,  the  expression  for  U  „  can  be  approximated 

—  w  O  \ 

by 


C4.ll) 


U  K  =  E  K  •  ’a (E  kE  ♦  EE  ^)  =  E  K  +  0(Eo)  ,  as  Eq  0  . 


(4.12) 


We  observe  that  with  U  and  U  K  given  by  (4.11)^  and  (4.12),  it  can  be 
verified  that  the  expression  (UU)  ^ =  2E  ^  and  thus  (4.11)2  and  (4.12) 
satisfy  (2.8)2  to  the  order  eQ.  Introducing  the  approximations  (4.11) 
and  (4.12)  on  the  right-hand  side  of  (3.15),  we  obtain 


*k  *  -  *k  -  !l>  -(!l®  !k5?,l 


+  W  ,,E  -  EE  „  -  2 [EE  (e„  0  e.)  -  (e.  0  eJE  ,E] 

j  N-w  n-wjN  *«•.  j  L  v  L  —  L/  ~  N  —  f  m-w 

+  2[E(eL  0  eK)E  L  -  E  L(fK  ®  eL)E]}  , 


(4.13) 


where  terms  of  0(eQ)  have  been  neglected. 
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4. 1  Estimate  for  the  angle  of  rotation 

With  the  help  of  (1.3),  we  can  estimate  the  norm  of  Aj,  from  the 
expression  (4.13): 

*  ||EE  K||  *2  II!!jl(!k®SL,II  *2  IIc!l®Sk!5,L?'I 

*  2||E(eLeeK)E(L|l  *2||EiL(!|teeL)E!|} 

<  6e1  +  13£q£,  *  ^(6  +  13£Q) 

Hence,  given  eo,  it  follows  from  (4.14)  that 

AK  =  0(£l)  as  e:  -  0  ,  (K  =  1,2,3)  . 

Prom  (3.14),  which  can  be  written  as  R  =  RA„,  the  norm  of  R  is  given  by 

•w  I  K  ‘'"•N  )1\ 

||R;KH2  =  tr{(R>K)TR>K}  »  ||  AkS|  2  . 

By  (4.16)  and  (4.14)  we  also  have  ||  R  J|  <£.(6  +  156  )  and  it.  follows  that 
the  order  of  magnitude  of  R  is  given  by 

**  t  ^ 

R>K  =  0(ej)  as  £x  -►  o  ,  (K  -  1,2,3)  . 

In  terms  of  (4.7)  with  m  =  3  and  n»l,  the  estimate  ||  R  „(|  noted  above 

-  >  »■ 

implies  (4.7)  with  k  =  2.  =  l,  and  we  conclude  that  quantities  of  0(6^)  are 
also  of  O(e^)  as  +  It  should  be  noted  that  although  (3.14),  which 
involves  R,  is  used  to  estimate  the  order  u f  magnitude  of  R  the 
result  (4.17)  is  independent  of  the  measure  of  smallness  £2* 


-  -  tj  ■,r  a1  timn't  \  r,  t~”  rrr: 


kTJisJSk 


usaii&Ex?j?3aitta 


(4.14) 


(4.15) 


(4.16) 


(4.17) 
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In  order  to  assess  the  effect  of  E  „  on  the  order  of  magnitude  of 

***  »  ^ 

the  angle  of  rotation,  we  now  use  the  representation  (2.13)  to  calculate 
the  left-hand  side  of  (3.14): 


RTR  k  =  (1  -  cos  6)  (y  ®  y  K  -  y  jK  ®  U) 

+  sin  9(1  -  cos  0)  { (y  K><  y)  ®  y  -y  ®  (y>K  *  y)> 

-  6  „e.  ®  (y  x  e. )  -  sin  0  cos  0  et  ®  (y  K  x  e. )  .  (4.18) 

,K~L  ~  ~J  i  ***-• 

The  above  expression  involves  both  the  angle  0  and  its  gradients.  In  order 
to  estimate  the  order  of  magnitude  of  the  latter  quantities,  we  proceed  to 
isolate  the  term  containing  9  „  in  (4.18) .  To  this  end,  we  first  observe 

i  N 

that  if  the  directions  y  „  vanish,  every  term  in  (4.18)  is  zero  except  the 

~  ,  K. 

one  containing  0  namely  -0  ^e^®  (y  x  e^)  •  Next,  let  ctj,  (K*  1,2,3)  be  a 

set  of  unit  vectors  along  the  directions  of  y  ^  t  0  and  choose  a  second  set 

of  unit  vectors  yv  (K=  1,2,3)  such  that  for  each  K  the  set  of  unit  vectors 

★  . 

(y,Yt,,cO  form  a  right-handed  orthogonal  triad  .  Then,  by  considering  the 

•w  —  -wK 

rr,  ★★ 

scalar  aR  •  RR  kyk  (no  sum  on  K) ,  it  can  be  verified  that 


*Our  preference  for  the  order  of  the  unit  vectors  in  the  right-handed  triad 
Oi»lK>ai()  is  siinP1>r  because  this  leads  to  a  positive  sign  on  the  right-hand 
side  o?  (4.19).  Alternatively,  a  negative  sign  would  result  on  the  right- 
hand  side  of  (4.19)  if  we  had  chosen  the  set  (P.oik'Ik)  as  a  right-handed 
orthogonal  triad. 

Details  are  given  following  (BIO)  in  Appendix  B. 
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5T?,kIk 


*-e,K(S“&*'Xi:)  -*6,k 


(no  Siam  on  K) 


From  comparison  of  (3.14),  (4.15)  and  (4.19)  we  may  now  conclude  that 
throughout  the  body  3  the  gradients  of  the  angle  of  rotation  have  the 
order  of  magnitude  given  by 


(4.19) 


0  K  =  Ofcj)  as  e:  ♦  0  (K-  1,2,3) 


(4.20) 


To  continue  the  discussion,  suppose  that  R =  I  at  some  material  point 
oX  e  3,  or  that  equivalently  0  =  0  at  X.  Then,  at  any  other  material 
point  X,  the  angle  of  rotation  can  be  calculated  from 


0 


rx 

~  O0/8Xk)Xk  dS 


where  dS  is  the  arc  length  of  an  arbitrary  curve  C  in  the  reference  con¬ 
figuration  of  3,  X  =  X^e^  is  the  unit  tangent  vector  to  C  and  the  integra 

* 

tion  is  performed  along  the  curve  joining  X  and  X.  By  means  of  the 

0~  -V 

usual  inequalities  for  integrals,  from  (4.21)  we  obtain  an  estimate 
for  the  magnitude  cf  0  in  the  form 


(4.21) 


lei  < 


f  X 

■  X 
o~ 


|39/3Xk| 


where  (0  ,)  and  X  are  the  values  of  0  and  X„  at  some  point  on  the  curve 

i  K  K  ,  K  K 

C,  respectively,  and  where  L  is  the  length  of  the  curve  C  joining  QX  and  X. 
Since  xT  are  components  of  a  unit  vector  and  since  the  estimate  (4.20) 

N 

holds  throughout  the  body,  from  (4.22)  we  conclude  that  the  angle  9  at 
every  point  of  3  has  the  order  of  magnitude 


(4.22) 


9  =  O(e^)  as  0  ,  (4.23) 

provided  that  R =  I  at  some  point  QX  e  3.  Recalling  the  power  series 
* 

The  existence  of  such  a  curve  in  3  is  ensured  by  the  connectivity  assump¬ 
tion  of  the  region in  section  2. 
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At- 


expansions 


a*'  a  “ 

sin  9  a  0  -  jp  +  •  •  •  ,  cos  0  =  l-  yj-+|y---*  > 

it  follows  from  the  estimate  (4.23)  that 

2 

sin  0  =  0(ep  ,  cos  9  =  1  +  0(ep  as  •+  0 

Using  (4.25)  in  the  representation  (2.13),  we  conclude  that  provided  R=  I 
at  X,  the  rotation  tensor  R  satisfies 

R  =  I  +  0(e1)  as  ■+  0  . 

Since  by  (4.5)  the  rotation  was  defined  as  RsI  +  OCe,),  it  fo1  lows  from 
(4.26)  that  provided  R  is  a  unit  tensor  at  some  point  qX  r.f  the  body, 
quantities  of  0(e2)  are  comparable  to  O(e^)  as  £^-+0. 

Thus  far  in  the  development  of  this  section,  the  order  of  magnitude 
of  E  „  (and  hence  also  of  R  or  0)  has  been  regarded  as  independent  of 

N  )  1\  V 

that  of  the  strain  E.  In  order  to  relate  the  estimate  of  the  rotation 
tensor  R  in  (4.26)  or  the  angle  of  rotation  0  in  (4.23)  to  the  order  of 
magnitude  eQ,  an  additional  assumption  must  be  made  concerning  the 
relationship  between  the  orders  of  magnitude  of  E  and  E  v  or  equivalently 
between  eq  and  e^.  But,  prior  to  such  an  undertaking,  we  need  to  dispose 
of  some  geometrical  preliminaries  and  definitions. 

Consider  now  any  unit  vector  v  and  rotate  this  by  the  rotation  tensor 
R  through  an  angle  a  resulti  ig  in  the  vector  R  v.  Since  the  proper 
orthogonal  tensor  R  is  length  preserving,  the  magnitude  of  v  remains 
unchanged  upon  rotation  and  we  have 

cos  a  =  v  •  R  v  , 

which  represents  the  projection  of  the  rotated  vector  R  v  along  v  (see  the 


(4.24) 


(4.25) 


(4.26) 


(4.27) 
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sketch  in  Fig.  1).  With  the  use  of  (2.13),  it  can  be  shown  from  (4.27) 

* 

that  cos  a  is  bounded  from  below  by  cos  9  in  (2.14),  i.e.  , 

cos  a  >  cos  9 

It  may  be  noted  here  that  the  angle  a  *  0  if  v  is  parallel  to  the  unit 

vector  y  in  (2.14),  while  a  *  0  if  v  is  perpendicular  to  y. 

To  continue  the  discussion,  let  8  stand  for  a  vector  defined  by  (see 

**• 

also  the  sketch  in  Fig.  1) 

8  =  R  v  -  v  or  R  v  =  v  +  8  . 

Keeping  in  mind  that  both  v  and  R  v  are  unit  vectors,  from  the  inner 
product  of  (4.29)  with  itself  we  arrive  at 

2 

8  •  8  =  Js|e|  =  1  -  v  •  R  v  =  1-  cos  a  , 

~  * W  A* 

2 

which  shows  that  k\&\  is  the  difference  between  unity  and  the  projection 
of  R  v  along  v.  Now,  let  <|>  denote  the  angle  that  the  vector  8  makes  with 
v  as  indicated  in  Fig.  1.  Then,  by  (4.29)^,  we  have 

8  •  v  =  1 8 1  cos  <j>  *  -  (1  -  v  •  R  v) 

Since  the  right-hand  side  of  (4.31)  is  -%|b|  by  (4.30)2*  it  follows  that 
cos  <j>  in  (4.31)  is  not  an  independent  quantity  and  is,  in  fact,  given  by 
cos  4>  =  -  H\&  \ . 

4 . 2  The  special  case  of  infinitesimal  kinematics 

Before  proceeding  further,  it  is  instructive  to  consider  the  case  of 
classical  infinitesimal  kinematics  in  which  both  the  strain  and  rotation 
are  small.  Thus,  we  introduce  the  following 
Details  of  the  argument  are  given  following  (Bll)  in  Appendix  B. 


(4.28) 


(4.29) 


(4.30) 


(4.31) 
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Definition  4.1.  Given  E=0(e_),  a  proper  orthogonal  tensor  R  is 

1  0  <v 

said  to  be  an  infinitesimal  rotation  with  respect  to  eq  if  for  any  unit 
vector  v,  the  vector  g  defined  in  (4.29)  satisfies 

3  =  0(0  as  e  -*■  0  . 

It  follows  at  once  from  (4.32)  and  (4.30)2  that  R  is  an  infinitesimal 
rotation  in  the  sense  of  (4.32)  if  and  only  if 

2 

1-vRv  =  0(e  )  as  £  +  0  , 

^  ^  0  0 

2 

i.e.,  the  projection  of  R  v  along  v  differs  from  unity  byO(£  )  as  £  -*■  0 
if  and  only  if  R  is  an  infinitesimal  rotation.  Clearly,  (4.33)  can  be 
used  to  state  an  alternative  definition  of  infinitesimal  rotation  which 
is  equivalent  to  (4.32). 

2 

Observing  from  (4 . 30) 3  and  (4.33)  that  cosa=l+0(£o)  as  eQ  +  0,  the 
inequality  (4.28)  together  with  the  fact  that  a  =  9  for  some  v  imply  that 

2  2  2 

cos  0  =  1  +  0(eq)  ,  sin  0  =>  0(eq)  ,  sin  0  ■  0(eq)  , 

0  =  0 (£  )  as  e+0 
v  o'  0 

Introduction  of  the  approximation  (4.34)  into  (2.14)^  ^  yields 

$  =  <t>T  -  0(e2)  ,  4,  =  n=  -QT=0(e)  as  £  -*■  0  . 

~  ~  o  ~  ~  0  0 

If  R  is  an  infinitesimal  rotation  with  respect  to  eQ  everywhere  in  the 
body,  then  in  terms  of  the  condition  (4.7)  we  have  m  =  2,  n  =  0,  k  =  £=l, 
and  quantities  of  0^)  are  also  of  0(eq)  as  e  -*-0. 


(4.32) 


(4.33) 


(4.34) 


(4.35) 
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It  should  be  emphasized  that  implicit  in  the  Definition  4.1  is  the 
assumption  that  the  order  of  magnitude  of  R  can  be  estimated  in  terms  of 
the  order  of  magnitude  of  E.  To  avoid  undue  complications,  we  have  post¬ 
poned  making  explicit  any  relationship  between  and  eq  but  the  conditions 
under  which  R  can  be  estimated  in  terms  of  E  will  be  examined  later.  In 
this  connection,  it  may  be  recalled  that  in  the  usual  kinematics  of  the 
infinitesimal  theory  the  restriction  of  smallness  is  imposed  on  the 
relative  displacement  gradient  H  defined  by  (2.6);  and  then,  it  follows 
at  once  that  both  the  strain  and  the  rotation  are  infinitesimal.  The 
approach  in  the  Definition  4.1  differs  from  the  usual  in  that  so  far  no 
restriction  has  been  placed  on  the  relative  displacement  gradient  H,  but 
infinitesimal  rotation  is  defined  with  respect  to  e  . 

In  the  context  of  infinitesimal  kinematics,  we  now  state  the  following 
Theorem  4.1.  The  relative  deformation  gradient  H  is  of  0(e  )  as 

e  -*-0  if  and  onlv  ifE=»U-I  =  0(e  )  ase  ->-0  and  R  is  infinitesimal 
0  '  ^  ^  ~  0  0  - 

rotation  with  respect  to  in  the  sense  of  Definition  4.1. 

Proof.  Let  the  tensor  U  be  specified  by  (4.11)2  and  suppose  that 
are  given  by  (4.35)  so  that,  after  the  neglect  of  terms  of  O(e^), 
the  rotation  R=I  +  0(eq)  as  e  ■*■0.  Then,  by  (2.1)  the  deformation  gradient 
F  =  I  +  0  (e  )  as  e  -*■  0  and  H»0(e  )  as  e  -*■  0 . 

Conversely,  if  H  defined  by  (2.6)  is  of  0(eQ)  as  eo  +  0,  then  from 
(2.3)2  and  (2.6) ^  we  have 
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(4 . 36) 


U2  -  I  +  H  +  HT  +  0(e2)  as  e  -*■  0  . 

A*  >y  >w  0  O 

2 

If  terms  of  0(eo)  are  neglected  as  e^O  in  (4.36),  we  have 
U  »  I  +  %(H+HT)  ■  I  +  0(e  )  as  e  ■»  0  . 

•w  «W  «W«mO  0 

Recall  now  that  R  »  Flf*  and,  with  the  help  of  (4.37),  obtain 

-w  «*** 

R-  F[I-^(H+HT)]  «  (I+H)  (I  -  -  ^HT) 

-  I  +  ?s(H-HT)  -  I  +  0(e  )  as  £  +  0  , 

•W  •«  X  0  0 

2 

where  terms  of  0(eQ)  have  been  neglected  as  e  -*0.  This  completes  th 3 
proof. 

It  may  be  worth  recalling  that  the  strain  gradients  E  ,,"0(e.)  as 

■V  J  h  <w  X 

•+■  0 ,  in  view  of  (4.2);  and  that  subsequently  the  use  of  this  measure  of 
smallness,  along  with  (4.18)  enabled  us  to  establish  the  estimate  (4.23) 
for  the  angle  of  rotation  8  provided  that  R  =  I  at  some  point  qX  £3.  The 
estimate  (4.23)  can  be  brought  into  correspondence  with  (4.34)^,  which 
results  from  Definition  4.1,  if  and  only  if  functions  of  O(e^)  are  assumed 
to  be  comparable  to  0(e).  In  terms  of  the  conditions  stated  immediately 
after  (4.4),  we  may  choose  k  =  1 ,  l  =  1  so  that  any  function  of 
O(e^)  =  0(tQ)  as  eo  +  0.  With  this  additional  assumption,  we  may  now  write 

E  v  «  0(e  )  as  e  •+•  0 

v  j  K.  O  0 

and  ttyis  is  consistent  with 

9  *  0(eo)  as  cq  -*•  0  , 

provided  R«I  a1  some  point  QX  €  3. 

We  are  now  in  a  position  to  state  the  following 


(4.37) 


(4.38) 


(4.39) 


(4.40) 
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Theorem  4.2.  Given  E  =  0(e  )  and  E  „>0(e  )  as  £  -*-0,  the  tensor  R 

1 mmmt 1  ^ **  —  —  0  v ;  *v  0  0  — 

associated  with  the  deformation  function  x  is  an  infinitesimal  rotation 
with  respect  to  £q  to  within  a  rigid  body  rotation. 

Proof.  Provided  that  R*  I  at  some  point  X  S  S>,  the  conclusion 
(4.40)  implies  that  R  is  an  infinitesimal  rotation  with  respect  to  e  . 

Thus  with  the  use  of  (2.13)  and  (4.35),  to  the  order  of>  approximation 
considered,  the  infinitesimal  rotation  is  given  by  (I  +  ft)  as  £  -*-0. 
Moreover,  by  Theorem  3.1  any  other  solution  of  (3.14)  differs  from  (I  +  ft) 
by  a  proper  orthogonal  tensor  function  of  time  (say  R  )  corresponding  to 
a  rigid  body  rotation  and  we  have 

R  »  R  (I +  ft)  as  e+0  , 

subject  to  the  condition  that  (I  + ft)  *  I  at  X  ■  or  equivalently 

ft  a  0  ,  R  *  R  at  X  =»  X  , 

■v  ..  ~  ~o  ~  o~ 

where  ft  =  0(:o)  as  £q-M)  is  defined  in  (4.35).  This  completes  the  proof. 

It  may  bo  observed  that  in  the  present  development,  apart  from  rigid 
body  rotation,  the  conclusion  ft  =  R  -  I  =  0(eQ)  as  eQ  +  0  is  a  derived  result. 
This  is  in  contrast  to  the  usual  approach  to  infinitesimal  kinematics, 
where  the  infinitesimal  nature  of  ft  (and  also  of  E)  is  implied  through  an 
assumption  of  smallness  imposed  on  the  displacement  gradient  H. 

The  restrictions  in  the  statement  of  Theorem  4.2  are  imposed  on  E  and 
E  (instead  of  on  the  usual  displacement  gradients)  and  for  an  explicit 
calculation  of  the  rotation  tensor  we  need  to  return  to  (3.14).  Alterna¬ 
tively,  in  view  of  the  representation  (4.41),  it  will  suffice  to  calculate 
the  skew- symmetric  tensor  ft  and  then  R  is  determined.  To  see  this, 
introduce  (4.39)  into  (4.13)  and  after  neglect  ot  terms  of  0(£p  as 

e  -*-0,  obtain 
o 


(4.41) 


(4.42) 
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(4.43) 


-K 


-AT 

~K 


E  ,(e„®  e.)  -  (e  ®  eJE  .  -  0(0 

*.,L.  N  ~  L  -<  L  -»  S  ,  L,  n.  U 


as 


£4-0 


But  the  expression  for  the  gradient  of  R  by  (3.14)  is  R  „  ■  R  A„  and 
substitution  from  (4.41)  results  in 


RJ1  \f  *  R  (I  +  * 

•'O  —  ***  "4  K 

Since  R  is  nonsingular  and  since  Q*0(e  )  and  A,,«0(e  )  as  e  ■+  0  by 
*.0  ~  o  ~K  -  o  o  4 

* 

(4.35),  and  (4.43),  respectively,  (4.44)  gives 

~,K  *  ~K  ' 

2 

where  terms  of  0(e  )  have  been  neglected.  The  tensor  ft  can  now  be 

•w  0 

determined  as  a  solution  of  (4.45)  subject  to  the  condition  (4.42). 

With  the  use  of  (2.1),  (4 . 1 1 ) 2  ani*  (4*41)  we  note  that  the  deforma¬ 
tion  gradient  is  now  given  by 


F  *  R  (I+ft  +  E)  as  e  *  Q  , 

—  *v  0  -V  «w  0 

2 

where  terms  of  0(gq)  have  been  neglected.  In  addition,  by  (2.6)  and 
(4.46)  we  also  have 


H  ■  (R„  -  I)  +  R  (ft  +  E)  as  e  4  o  . 

•W  —0  M  -wO  <W  •»  0 

In  the  above  expression,  it  is  H  -  (R  -  I)  and  not  H  itself  which  is 
infinitesimal  with  respect  to  eo> 

4. 3  Small  strain  accompanied  by  moderate  rotation 

The  foregoing  development  [between  (4.32)  and  (4.47)]  which  began 
with  Definition  4.1  dealt  with  infinitesimal  kinematics.  We  now  return 
to  our  main  objective  and  introduce 

Definition  4,2.  G  ven  E  =  0(eQ),  a  proper  orthogonal  tensor  R  is 
*  _  ’  ' 

The  component  form  of  Eqs.  (4.45)  are,  of  course,  the  same  as  those  used 
in  the  infinitesimal  theory  of  notion  mentioned  in  the  first  paragraph 
of  section  1. 


(4.44) 


(4.45) 


(4.46) 


(4.47) 


31 . 


1 


said  to  be  a  moderate  rotation  with  respect  to  e  if  for  any  unit  vector 

°  ] 

v,  the  vector  3  defined  in  (4.291  satisfies  ? 

~  •  | 

k  i 

8  *  0(0  as  e  *  0  .  (4.48)  i 

-^00  ] 

It  follows  at  once  from  (4.48)  and  (4.30)2  that  R  is  a  moderate  rotation  j 

in  the  sense  of  (4.48)  if  and  only  if  [compare  with  (4.33)]  1 


l 


l ' 

< 

j ‘ 

i 

I 

i'' 


i 


i 


t* 


1  -  v  •  R  v  ■  0(eQ)  as  £q  ■+■  0  ,  (4.49) 

i.e.,  the  projection  of  R  v  along  v  differs  from  unity  by  0(e)  as  £+0 
in  the  case  of  moderate  rotation. 

Observing  from  (4 . 30) ^  and  (4.48)  that  for  moderate  rotation: 
cos  a  ■  1  +  0(e  )  as  eo+0,  the  inequality  (4.28)  together  with  the  fact 
that  a  a  0  for  some  v  imply  that 

cos  6  *  1  +  0(eQ)  ,  sin"  0  a  0 (Eq)  >  sin  9a0(e^)  ,  9“0(e^)  as  £q  ■+•  0  .  (4.50) 

Clearly,  (4.49)  can  be  used  to  state  an  alternative  definition  of 
moderate  rotation  which  is  equivalent  to  (4.48). 

As  in  the  Definition  4.1,  wo  again  observe  that  implicit  in  the 
Definition  4.2  is  the  assumption  that  the  order  of  magnitude  of  R  can  be 
estimated  in  terms  of  E.  Again  to  avoid  undue  complications,  we  have 
postponed  making  explicit  any  relationship  between  and  eQ  but  the 
conditions  under  which  moderate  rotation  can  be  estimated  in  terms  of 


i 


L 


infinitesimal  strain  will  be  examined  later.  However,  two  aspects  of 
the  conclusions  (4.50)  for  moderate  rotation  may  be  noted  here: 

(1)  The  results  (4.50)  have  been  obtained  with  respect  to  and  are 
independent  of  defined  by  (4.2);  and  (2)  the  angle  of  rotation  in 
(4.S0)  is  of  0(ep  in  constrast  to  that  of  0(eQ)  in  (4.34)  for  infinitesi¬ 
mal  rotation.  Thus,  after  substituting  from  (4,50) ^  ^  in  (2.14)^  ^  and 


3/2 

neglecting  terms  of  0(eQ  ),  the  functions  T  and  $  in  (2.13)  for  moderate 
rotation  are  approximated  according  to 


T  «  -  4/T  -  0(A  ,  4>  *  $T  *  0(e  ) 

—  —  O  *«-  —  0 

T  T  ■  2$  *  0(e  )  as  e  -*■  0 
***  0  0 


(4.51) 


It  should  be  noted  here  that  if  R  is  a  moderate  rotation  with  respect  to  e 

-  o 

everywhere  in  the  body,  then  in  terms  of  the  condition  (4.7)  with  m»2,  n»0, 
k*2,  £-1,  quantities  of  0(e_)  are  also  of  0(e"4)  as  e  +0. 

L  0  0 

We  now  state  the  following 

Theorem  4,3.  IfE»U-I*0(e  )  as  e  -*-0  and  R  is  a  moderate  rotation 

A*  0  0  *V 

in  the  sense  of  Definition  4.2,  then  the  tensor  H  defined  by  (2.6)  is  of 
0(e*)  as  e  +  0.  Conversely,  if  H»0(e4)  as  e  +0  and,  in  addition,  if 
the  symmetric  part  of  (2.6),  i.e.,  H  +  is  of  0(eq)  as  £q -*■(),  then 

U-I»0(eo)  and  R  is  a  moderate  rotation  with  respect  to  eQ. 

Proof .  We  first  prove  the  first  part  of  the  theorem.  If 
U  -  I*E"0(eo)  as  eo  +  0  and  R  is  a  moderate  rotation  given  by  (2.13) 
with  <t>  and  V  estimated  by  (4.51),  then  recalling  (2.1)  we  have 


F  a  I  +  f  +  <J>  +  E  ,  E  ■  0(e  )  , 

M  m*  «■>  <V  N  0 

given  by  (4.51)  as  e  +  0 

and  H  >  0(e|j)  by  (2.6) . 

We  now  turn  to  the  converse  part  of  the  theorem.  By  assumption, 

\f  T 

H»0(e'J)  and  H  +  H  ■  0(eq)  as  eo  +  0.  Hence,  substitution  into  (2.7)  ^ 
yields  E  »  0(eQ)  as  e  -*■  0  and  by  (4 . 11) ,  we  also  have  U  «  I  +  0(eQ) . 
Further,  in  order  to  obtain  the  desired  estimate  for  R,  from  (2.1)  we 

N« 

have  R=  F  U  *  and  from  (4.11)2  we  recall  that 

U'1  -  I  -  E  -  I  -  Jfi(H  +  HT+  HTH)  as  e  •*>  0  . 


(4.52) 


(4.53) 


Then,  R  can  be  expressed  as 


•Uw.-mUj: iii iiii ^ ■  /. -  gb .... . a ...  ,  ...a  i.y,  - .  .  ,t..4  .. 


L*iL!t<  U//1AV  iL'l.iU'A  i'tl 


33. 


C4.54) 


R  *  (I  +  H) {I  -J*(H  +  HT  +  HTH)} 

=  I  +  J*(H  -  HT)  -J$HTH  +  0(e3/2)  . 

A  close  examination  of  (4.54)  reveals  that  if  we  identify  the  skew- 
T  T 

symmetric  J$(H  -  H  )  as  4*  and  the  symmetric  -  H  as  $  so  that  these 
quantities  can  meet  the  condition  (4.51)_,  then  R  in  (4.54)  can  be 

w 

identified  as  moderate  rotation.  To  show  this,  we  recall  that  by 
assumption 

H  *  -  HT  +  0  (e  )  and  HT  3  -  H  +  0  (e  )  as  e  -*-0  . 

•W  N  IS*  0  "■*  **  *■>  0  0 

Thus,  after  forming  the  product  !  !  »  [^(H  -  H^) ] [^(H  -  HT) ]  and  neglecting 
3/2 

terms  of  0(e  ) ,  we  obtain 

«v  O 

[Jj(H-HT)][Jj(H-HT)]  -  ?^{-2HT  +  0(e  ) }{2H  +  0(e)) 

T 

3  -  H  H  , 

which  meets  (4.51) ^ •  Hence,  R  given  by  (4.54)  is  a  moderate  rotation  and 
the  proof  of  the  theorem  is  complete. 

We  recall  once  more  that  the  strain  gradient*  E  ^O(e^)  as  e^  +  0, 
in  view  of  (4.2),  and  that  the  estimate  (4.23)  for  the  angle  of  rotation 
0  is  also  obtained  with  respect  to  e^.  Tins  estimate  (4.23)  can  be 
brought  into  correspondence  with  (4.50)^,  which  results  from  Definition 

4.2,  if  and  only  if  functions  of  Q(ep  are  assumed  to  be  comparable  to 

1 , 

O(e’) .  In  terms  of  the  conditions  stated  following  (4.4),  we  may  choose 

u 

k*2,  i  ■  1  so  that  any  function  of  O(e^)  »  O(e^)  as  e  +0.  With  this 
additional  assumption,  we  may  now  write 

E  =  0(c'3)  as  e  -*•  0 
~,K  ~  o'  o 


(4.55) 


(4.56) 


(4.57) 


34. 


and  this  is  consistent  with 


6  *  O(e^)  as  £0  -*•  0  , 


(4.58) 


provided  that  R  =  I  at  some  point  qX  £  SB. 

We  now  state  a  theorem  concerning  the  order  of  magnitude  of  R: 

k 

Theorem  4.4.  Given  E  =  0(e  )  and  Grad  E»0(e*)  as  £  -*-0,  the  tensor  R 

«v  o  o  0  *v 

associated  with  the  deformation  function  x  is  a  moderate  rotation  with 

respect  to  eQ  to  within  a  rigid  body  rotation. 

Proof.  Provided  that  R=I  at  some  point  the  conclusion  (4.58) 

implies  that  R  is  a  moderate  rotation  with  respect  to  £q  in  the  sense  of 

Definition  4.2.  Thus,  with  the  use  of  (2.13)  and  (4.51)  and  to  within 
3/2 

terms  of  0(eo'“),  the  moderate  rotation  is  given  by 

I+$+y  as  e  ■+•  0  , 

~  ~  ~  o  * 

4>  and  T  specified  by  (4.51) ^  2 

Moreover,  by  Theorem  3.1  any  other  solution  of  (3.14)  differs  from  (4.59) 
by  a  proper  orthogonal  tensor  function  of  time  (say  Rq)  corresponding  to 
a  rigid  body  rotation  and  we  have 


(4.59) 


R  =  Rq(I  +  *  +  T)  , 


$  and  T  specified  by  (4.51).  _ 

-  -  L  >  / 


subject  to  the  conditions  that 


(4.60) 


*  -  V  -  0  ,  R  =  R  at  X  =  X 
-  ~  -  -  ~o  ~  o~ 


(4.61) 


and  the  theorem  is  proved. 

It  should  be  emphasized  here  that  the  approximations  for  moderate 
rotation  occur  only  in  the  functions  $  and  T  as  specified  by  (4.51)^  ^ 


The  rotation  R  in  (4.60)  may  in  fact  be  large,  in  view  of  the  presence  of 

Rq  which  represents  a  rigid  body  rotation.  Indeed,  it  is  the  proper 

T 

orthogonal  quantity  R  R  which  may  be  termed  moderate  rotation  and  not  R. 

•wO'V  *v 

The  restrictions  in  the  statement  of  Theorem  4.4  are  imposed  on  E  and 
F.  v  and  for  an  explicit  calculation  of  the  rotation  tensor  we  need  to 
return  to  (3.14).  Alternatively,  in  view  of  the  representation  (4.60),  we 
may  calculate  $  and  4*  to  the  order  of  approximation  considered  and  then  R 

is  determined.  To  show  this,  introduce  (4.57)  into  (4.13)  and  after  neglect 
3/2 

of  terms  of  0(e  )  as  e  +0,  obtain 

—  0  0 


i 


as  eo"°  ■ 

Substitution  of  (4.60)  into  R  =R  \v,  which  is  obtained  from  (3.14), 

"*  )  IN  i\ 

results  in 


(4.62) 


5o», *  WfPiic  • 


But,  with  the  use  of  the  order  of  magnitude  estimates  (4.51).  7  and  (4.62), 

1 

3/2 

after  neglecting  terms  of  0(eq  )  and  remembering  that  Rq  is  nonsingular, 
from  (4.63)  we  obtain 


(4.63) 


U*Um  Ik'Hk 


(4.64) 


Since  'V  v  represents  the  skew- symmetric  part  of  (4.64),  we  may  write 
~  >  K 


T  .  «  ki (<*>  „  +  T  „)  -  (4>  .  +  «P  J  ) 

**  )  A  •<  j  N  *4  j  N  j  N 


=  (^K  + 1  tt?  ’  (~K  + 1  Ak)  > 


=  Ak  +  WA|(-A1(?)  , 


(4.65) 


where  (4.62) ^  has  been  used  in  obtaining  (4.65) Given  E  ^  (and  hence  A^) , 
(4.65) ^  may  be  viewed  as  the  differential  equation  for  T ;  and,  once  ¥  is 


. u-  f 


determined,  4>  can  then  be  found  from  (4.51) y 

It  is  useful  to  record  here  the  expressions  for  the  deformation 
gradient  F  and  the  relative  deformation  gradient  H  associated  with  moderate 
rotation.  Thus,  with  the  use  of  (2.1),  (4 . 11) 2 ,  (4.60)  and  upon  neglect  of 
terms  of  0(e^2),  the  deformation  gradient  can  be  expressed  as 

F  =  R  (I  +  $.+  ¥  +  E)  as  E  0  . 

_  ~o  ~  ~  ~  ~  o 

Similarly,  by  (2.6)  and  (4.66),  we  have 

H  =  (R0  ‘  P  +  So  +  J  +  P  ■ 

In  view  of  the  remarks  made  in  the  paragraph  following  (4.61),  it  should 
be  noted  that  in  the  above  expression  it  is  [H  -  (R  -  I)]  and  not  H  itself 
which  is  moderately  large. 

4.4  Construction  of  a  motion  which  results  in  a  properly  invariant  small 
strain  and  moderate  rotation. 

From  among  all  particles  of  S,  let  Y  be  chosen  as  a  pivot  whose  trans¬ 
lation  and  rotation  is  specified  (Casey  and  Naghdi  1981).  Then,  cor- 

★ 

responding  to  any  motion  x>  we  can  construct  another  motion  X  by  removing 
from  x  the  translation  and  rotation  of  the  pivot.  Thus,  we  can  write 

x*  =  X*(X,t*)  =  R^X(X.t)  -XLX.t)}  +  X  , 

•s.  —  —  *»o  —  -<*  ~ 

★ 

t  =  t  -  c  , 

where  qX  denotes  the  position  vector  of  the  pivot  Y  in  the  reference  con¬ 
figuration  <  ,  Ro  =  R(QX,t)  is  only  a  function  of  time  and  c  is  a  real 

★  ★  * 

constant.  The  configuration  of  &  at  time  t  in  the  motion  X  is  <  .  We 
observe  that  (4.68)  is  of  the  form  (2.9)  with 

Q(t)  =  Rq  ,  a(t)  =-  Rjx(0X.t)  +0X  ,  a  =  -  c  . 

Then,  with  the  use  of  (2.10),  (2.11),  (4.60)  and  (4.66),  we  easily  conclude 
that 


(4.66) 


(4.67) 


(4.68) 


(4.69) 


37. 


F  sI  +  Q  +  ^  +  E  as  £  -*■  0  , 

~  ~  ~  ~  o 


R=I  +  $  +  4'  as  e->-0  , 

-  ~  -  o 


and  that 


(4.70) 


C=C,U=U,  E=E,E=E.,. 

~  ~  K  -,K 


(4.71) 


In  the  configuration  k  ,  the  relative  displacement  u  =X  -X,  and  the 


*  T 

rotation  tensor  R  =  R  R  are  such  that 
-  ~o~ 


U*(  X,t*)  =  0  ,  R* (  X,t*)  =  I  . 

~  u~  —  ~  o—  *** 


(4.72) 


The  expressions  for  Grad  u  =  u  „  ®  ev  are  again  of  the  form  (2.6)  but  with 

*»  j  1\ 

H  replaced  by 


H  =  u  ®e=$  +  'K  +  E  . 


(4.73) 


It  is  clear  from  (4.51).,  (4.11)  and  the  right-hand  side  ot  (4.73),,  that 

*  k  * 

u  .=0(0  as  e  -*-0.  From  this  result  and  the  fact  that  u  vanishes  at 
~,K  o'  o 

qX  by  (4 . 72) ^ ,  it  follows  that  the  displacement 


u  (X,t  )  =  0(ep  as  eQ  -*■  0 


(4.74) 


throughout  the  body  ,  Similar  to  (2.7)^,  the  strain  measure  E  can  be 


calculated  in  terms  of  H  so  that 


★  *  ★rr 

c,  =  %(H  +H  +H  H) 


■  *  Sk  ®  s!k *  (e*K  •  s!l>!k  ®  El1  • 


(4.75) 


It  should  bp  noted  that  in  the  above  formulae  for  infinitesimal  strain. 


The  line  of  argument  here  is  similar  to  that  employed  between  (4.20)- 
(4.23). 


38. 


44*"  -  '.I  .  — ~  r j. ; -I:, .  '  7.:,.',;’: .  -liy 


it  ^r^ri  **rt  -4  it  y 

while  both  quantities  H  +H  and  H  are  of  0(e  ),  H  =  0(e4)  as  e  -*-0. 

~  ~  ~  -  ~  0  ~  ~  0  0 

It  has  been  shown  in  the  paper  of  Casey  and  Naghdi  (1981,  Theorem  3.1) 
that  two  motions  ,x  and  _x  o£  &  differ  by  a  rigid  motion  if  and  only  if 

1**  Z  — 

*  *  ★ 

iX  = 0X  »  i-e.,  by  construction  x  remains  unchanged  if  x  is  replaced  by 

4,  4.  >r 

X  in  (2.9).  Thus,  with  (x  )  defined  by  [see  section  3.2  of  Casey  and 
Naghdi  (1981)] 

(xV(X,t)  =  OO^VcX’t*)  -X+(X,t+)}  +  X  (4.76) 

and  with  c  chosen  equal  to  -a  in  (2.9),  all  quantities  with  respect  to  the 
★ 

configuration  k  remain  unaffected  by  an  arbitrary  (not  necessarily  small) 
superposed  rigid  body  motions  and  we  conclude  that 

+  *  *  +*  * 

CX  )  =  X  ,  P  =  F  =  I  +  « *  l  +  E  , 

R+*  =  R*  ■  I  +  $  +  ¥  ,  H+*  =  H*  =  $  +  '}'  +  E  ,  (4.77) 

4.  ★  *  x*  4.*  *  4.*  * 

C  =  C  ,  U  =  U  ,  E  =  E  ,  Er«E  -  . 

~  ~  —  ~  —  —  «-  j  ^ ,  N 

+  *  +  * 
where  (F  )  =  Grad(x  )  . 


5 .  Invariance  of  constitutive  equations  with  small  strain  accompanied 
by  moderate  rotation 

We  briefly  discuss  here  the  invariance  of  constitutive  equations  (and 

hence  that  of  a  complete  theory)  for  small  strain  accompanied  by  moderate 

rotation.  Although  the  development  of  this  section  is  specifically 

carried  out  for  the  case  of  an  elastic  material,  it  will  be  clear  that 

our  main  conclusion  reached  holds  for  any  material.  The  procedure  for 

constructing  a  properly  invariant  theory  in  the  presence  of  moderate 

rotation  is  similar  to  that  used  by  Casey  and  Naghdi  (1981)  for  an 

infinitesimal  theory  of  motion  in  which,  apart  from  superposed  rigid 

body  motions,  both  strain  and  rotation  are  small. 

The  notations  for  the  mass  density,  in  the  configuration  k,  the 

outward  unit  normal  to  the  surface  the  stress  vector  t  acting  on 

3iR  and  the  Cauchy  stress  tensor  T  were  introduced  in  section  2.  We 

* 

denote  the  corresponding  quantities  in  the  motion  x  »  introduced  in 

*  ★  *  *  * 

(4.68),  and  in  the  configuration  <  by  p  ,n  ,t  and  T  ,  respectively. 

-w  ~  ** 

T  T 

Recalling  that  R  (QX,t)  =  Rq  in  (4.68)  is  a  function  of  time  only  and  plays 
the  role  of  Q(t)  in  (2.9),  as  noted  also  in  (4.69) it  follows  from 
(2 . 16) x ,  (2 . 16) ^ >  (2.17)  and  (2.18)  that 

p  =  o  ,  n  =  kTf  X,t)n  ,  t  =  r'(  X,t)t  , 

—  ~  0  **  ~  ~  ^  0>s>  ~ 

(5 

1*  =  5  • 

+  * 

Similarly,  associated  with  the  motion  (x  )  defined  by  (4.67),  we  have 
+  *  +  *  +*  +* 

the  quantities  p  ,  n  ,  t  and  T  .  These  quantities,  with  the  help  of 
(5.1).  0  ,  .,  (2.11)  and  (2. 16) -(2. 18)  transform  according  to  (Casey  and 
Naghdi  1981,  section  3): 


* 


an strain .« in-re":.-;".'™ 


YVk?A.t  . 


+*  *  +*  *  +*  * 
p  ■  p  ,  n  3  n  ,t  =t  , 

T+*  3  (Q(t)R  }TQ(t)T  QT(t){Q(t)R  )  =  T* 

—  «v  ^  ■%.  —  iv  vO  — 


(5.2) 


For  an  elastic  material,  let  \|)  =  ^(X,t)  denote  the  elastic  strain 

+  * 

energy  per  unit  mass  in  the  configuration  k,  Also,  let  \J>  and  if)  denote 

+  * 

the  strain  energy  per  unit  mass  in  the  configurations  k  and  <  ,  respec- 

+  * 

tively.  We  assume  that  ^  and  it  then  follows  that,  ip  =  4  The 
nonlinear  behavior  of  an  elastic  solid  may  be  characterized  by  the 
constitutive  equation 


(5.3) 


where  ■  ip (E) .  We  observe  that  in  view  of  (2.10),  (2.16)^  and  the  fact 

that  E+  3  E  the  value  of  T+  of  the  stress  tensor  given  by 

(5.3)  for  the  motion  X+  satisfies  (2.18),  so  that  (5.3)  is  a  properly 

•k 

invariant  constitutive  equation.  The  Cauchy  stress  T  in  the  motion 

Ilf 

X  has  the  form 

T* .  %py <4  •  (4)t>  (fV 

~  3E  3E 


and  satisfies  (5.2)^. 

i. 

Suppose  now  that  the  motion  is  such  that  F.  »0f£  )  and  E  „*--0(e'“)  as 
£q  +  0  and  recall  that  for  such  a  motion  the  deformation  gradient  is  given 
by  (4.66).  It  then  follows  from  the  local  equation  for  conservation  of 
mass,  namely  p  det  F*po,  that  for  small  strain  p  can  be  approximated  as 

p  =  pQ  +  0(eQ)  as  £o  -*•  0  ,  (5.5) 


where  pQ  is  the  mass  density  in  the  reference  configuration  kq.  Further, 
assume  for  simplicity  that  the  strain  energy  \p  is  quadratic  in  E  so  that 
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. ..... 


ir  r  r  'ir 


J-iiAd .  j  . 


It  is  now  clear  that  the  constitutive  equation  (5.8)  meets  the  require- 

t  ment  (5.2).  and  is  properly  invariant.  As  noted  in  subsection  4.4,  the 

[  ^  4 

f  method  of  construction  of  Casey  and  Naghdi  (1981) ,  which  is  also  used  in 

I 

j  the  development  of  this  section,  removes  from  all  motions  the  translation 

[■  |  -  and  rotation  at  any  particle  Y  of  the  body  called  a  pivot.  But,  even  in 

the  presence  of  moderate  rotation,  it  can  be  demonstrated  that  it  does  not 
matter  which  particle  is  chosen  as  pivot.  Indeed,  it  is  shown  in  Appendix  C 
that  the  theory  in  which  the  deformation  gradient  given  by  (4.66)  is  con- 

l  j 

structed  with  Y;  as  pivot  coincides,  to  within  terms  of  ),  with  that 

j  having  Y  as  pivot. 

■ 

r 


Appendix  A 


A  brief  account  of  notation  and  mathematics?,  terminology  used  in  the 
paper  is  given  at  the  end  of  section  1.  In  this  Appendix  we  collect 
additional  terminologies  and  mathematical  results  utilized  in  the 
development  of  the  paper. 

Any  linear  mapping  between  a  set  of  vectors  and  a  set  of  second  order 
tensors  will  be  regarded  as, a  third  order  tensor.  In  particular,  the 
tensor  product  a  ®  b  ®  c  of  any  three  vectors  a,b,c  £  V  is  a  third  order 

tensor  defined  by  (a  ®  b  &  c)v  *  c  •  v  a  ®  L,  and  we  also  define  a  product 

(a  ®  b  ®  c)  [u  ®  v]  *  (b  •  u) (c  •  v)a,  for  any  vector  u,v  e  V.  We  note  that 

the  tensor  product  between  a  second  order  tensor  and  a  vector,  i.e., 

T  ®  v  or  v  ®  T  for  any  second  order  tensor  T  and  any  vector  v,  is  also 

defined  as  a  third  order  tensor.  The  product  of  a  third  order  tensor  and 

a  second  order  tensor  is  again  a  third  order  tensor  defined  by 

(a  ®  b  ®  c)  (u  ®  v)  =  c*.  *  u(a  ®  b  ®  v)  =  (a  ®  c)  (u  ®  b  ®  v)  , 

(Al) 

(a  ®  b  ®  c)T  =  a  ®  b  ®  (TTc)  ,  T(a  ®  b  ®  c)  =  (Ta)  ®  b  ®  c  , 

for  any  vectors  a,b,c,u,v  6  V  and  any  second  order  tensor  T.  Since  the 
transpose  of  a  third  order  tensor  is  not  uniquely  defined,  we  introduce 
three  kinds  of  transpose  distinguished,  respectively,  by  superscript 

VVV  ™us 

(a  ®  b  ®  c)  =  a  ®  (b  ®  c)  =  a®c®b  , 

T 

(a  ®  b  ®  c)  ^  =  (a  ®  b)^  ®c  =  b®a®c  ,  (A2) 

T3 

(a&b®c)  =c®b®a 

for  any  vector  a,b,c  £  V,  and  we  note  that  the  definition  (A2)j  involves 
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the  transposition  of  c  and  a.  It  is  clear  that  as  a  consequence  of  these 

definitions,  our  definition  for  the  transpose  can  be  related  to  those 

for  the  transpose  T.  and  T-.  Keeping  this  in  mind,  it  is  convenient  to 

T3  ti  T,  T  T  T  T2 

introduce  the  abbreviation  <-A  =  [(._>!  )  ]  »  [0_4  )  J  for  a  third  order 

tensor  Jl. 

The  linear  napping  from  the  set  of  second  order  tensors  into  itself 

is  a  fourth  order  tensor.  In  particular,  the  tensor  product  a  &  b  8  c  ®  d 

of  any  four  vectors  a,b,c,d  £  V  is  a  fourth  order  tensor.  It  is  useful  to 

record  the  relationship  (a  &  b  c  ®  d)  [u  ®  v]  *  (c  •  u)  (d  •  v)a  ®  b,  which 

is  a  second  order  tensor.  The  transposed  of  a  fourth  order  tensor  >K  is 

T 

defined  by  the  relationship  B*JC[A]  *  A  *.X  [B]  for  all  second  order 
tensors  A,B.  Clearly,  (a  8  b  %  c  S  d)*1"  =  (c  ®  d  ®  a  <8  b)  . 

We  discussed  the  component  representation  of  second  order  tensors 
in  section  2.  Similarly,  for  a  third  order  tensor  Jl  or  a  fourth  order  tensor 

•5"e  “ f  “rite  d  3,111  •?  '  JtKLMNjK®SL®SN®?N'  "here 

JKLM  ■  SK‘  andJC|aMN  ‘  -  ^M®!^' 

The  gradient  of  a  scalar-valued  function  <j>(X)  may  be  written  as 


Grad  <J> 


(A3) 


Similarly,  the  gradient  of  a  vector-valued  function  v(X)  and  the  gradient 
of  a  second  order  tensor-valued  function  T(X)  will  be  denoted  by 
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3v. 

Crad  r.  ■  5^  'l 

^LM 

Grad  T  ■  tts —  e.  ®  e„  ® 
**  9X^  -L  ~M 

while  their  components  can  be  written  as 
(Grad  v)eM  ■  (v^  ®  ;K).;H 
(Grad  T)eH-  (T,«  ;K).M 
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Appendix  B 


This  appendix  provides  details  of  the  calculations  of  several  results 
and  formulae  stated  in  sections  3  and  4.  We  begin  by  noting  that  the  con¬ 
ditions  (3.10)  or  (3.11)  are  regarded  .  s  the  system  of  differential  equa¬ 
tions  for  the  determination  of  F  from  the  knowledge  of  C.  As  such  they 
were  obtained  as  necessary  conditions.  To  show  that  they  are  also 
sufficient,  we  consider  the  expression  of  the  form 


(FiLFiM},K  *  FiL,KFiM  +  FiLFiM,K  *  (B1) 

Assume  that  F,.  are  the  solutions  of  (3.10),.  Then,  by  direct  substitution 

XL  b 

of  (3.10)2  on  the  right-hand  side  of  (B1 ) ,  it  follows  at  once  that 


(FiLFiM},K  r  CLM,K 


(B2) 


Thus,  any  solution  F^  to  (3.10)  satisfying  the  condition  that  ^LFiMaClM 
at  one  point  will  have  the  property  that  *  CLM  evoryw,lere •  Moreover, 

since  the  right-hand  side  of  (3.10)  is  symmetric  with  respect  to  the  indices 
(K,M),  it  follows  that  a  solution  to  (3.10)  will  also  satisfy  (3.6)  and 
this,  in  turn,  ensures  the  existence  of  the  deformation  function 

s» 

It  was  remarked  following  (3.11)  that  since  U  can  be  determined 
uniquely  from  the  knowledge  of  C,  with  the  use  of  (2.1),  we  obtain  the 
differential  equations  (3.12)  for  the  determination  of  the  rotation  R. 

To  show  this,  we  substitute  (2.1)  into  the  expression  on  the  left-hand  side 
of  (3.10)  and  obtain 


FiLFiM,K  "  RiNUNL(RiPUPM) ,K  "  UNLUNM,K + UNLRiNRiP,KUPM  ’ 


or  after  rearrangement  of  terms 


RiMRiN,K  =  UML(FiLFi)',K  '  UQLUQP.K}UPN 


(B3) 
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Next,  we  substitute  from  (3.10)  on  the  right-hand  side  of  (B3) ,  make  use 

°£  CKL"UKMUML  by  (2-3)2  along  with  S^.k"  ^iM^N.K^M®  «N  Snd  °btain 
(3.14)  and  (3.15).  The  component  forms  of  these  equations,  namely 


!T!,K  -  Vk)UU 

*  ump^u(JP,l  *  'V,p")qkuln'Sm  ®  2n 


correspond  to  Eqs.  (8) -(10)  in  Shield's  paper  (1973). 

To  obtain  the  necessary  and  sufficient  condition  for  the  existence  of 
R,  we  first  recall  that  the  necessary  and  sufficient  condition  for  the 
existence  of  the  deformation  gradient  tensor  F  which  satisfies  (3.11)  is  that 


F,KL  “  F,LK  0r  FiM,KL  "  FiM,LK  * 

Again  substitute  from  (2.1),  and  after  cancellation  of  identical  terms, 
(B5)  becomes 


(B5) 


R  KTU  +  R  U 


~,LkH +  2  y,LK 


Clearly  for  a  given  U,  U  v.  ■  U  . v  and  since  U  is  nonsingular,  it  follows 

•v  ~  |  I\ L  j  LjK  «w 

from  (B6)  that  the  necessary  and  sufficient  condition  for  the  existence 
of  R  is 
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This  corresponds  to  Eq.  (11)  in  Shield's  paper  (1973),  which  was  obtained 
as  a  special  case  of  more  general  results  discussed  by  Thomas  (1934) . 

In  order  to  verify  the  truth  of  the  inequality  (4.4),  we  first 
observe  that  for  real  numbers  a  and  b,  0<a<b,  the  following  two 
inequalities  hold: 

1_  1_ 

an  <  bn  ,  an  <  bn  .  (B9) 


Next,  with  the  help  of  (B9)2j  from  the  result  (4.3)  which  holds  for  k,£, 
positive  integers,  we  obtain  (e^)^^<  (Ce^)*^  or  and  hence 


<  c^1/k£^/k3n 


(BIO) 


Then,  from  the  condition  on  II h.  (E  „)||  given  following  (4.2),  namely 

1  ~  )  N 

1 1 h -  ( E  „)  || <  De1?  as  e.,+0,  the  inequality  (4.4)  follows  from  (BIO). 

To  verify  the  truth  of  (4.19),  with  the  use  of  (4.18)  consider  the 
T 

scalar  ct  •  R  R  vyv  and  observe  the  identities  (no  sum  on  K) 

aK  •  (p  ®  M>k)Tk  b  V  *  *  Yk)  » 

ojK  •  t(M>Kx  y)  *  m]yk  *  u  •  YK(y  Kx  v  ■  aK)  » 

V!l®  =  '  Sl3  cy,K x  Sl  *  Ik5  *  H,kx  <V  Sl^l'Ik 

sPfxa  'y„  , 

V!i®  c>ix  Sl’Ik  *  (2k  ’  !l>  <H  *  ?l  '  Ik>  '  H  *  “k  ‘  Pc  • 

Since  (p.y^.a^)  form  a  right-handed  orthonormal  triad,  (B1 1) ^  2  vanish 

identically  and  (B1 1 ) ^  is  equal  to  -1.  In  addition,  since  is  directed 

T 

along  p  „ , (B 1 1 )  _  vanishes  also  and  the  scalar  a  •  R  R  vyv  reduces  to  (4.19) 

- , N  O  —  — ) N-N 

We  next  establish  the  validity  of  the  inequality  (4.28),  where  6  is 
the  angle  of  rotation  for  the  tensor  R,  and  a  is  the  angle  between  a  unit 


(Bll) 


hBdIEM£.-£iHiZ2l 
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vector  v  and  the  rotated  unit  vector  R  v.  With  the  help  of  the  representa¬ 
tion  (2.13),  we  have 

2 

cos  a  ■  v  •  (R  v)  =  (y  •  v)  (1  -  cos  8)  +  cos  9  .  (B12) 

2 

Observing  that  (1  -  cos  9)>0,  and  (y  •  v)  >0,  we  conclude  from  (B12)  that 
cos  a>cos  0  and  (4.28)  is  verified. 
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Appendix  C 


In  a  recent  paper,  Casey  and  Naghdi  (1981)  have  shown  that  a  theory 
in  which  both  the  deformation  and  deformation  gradient  are  small  can  be 
constructed  so  as  to  possess  desirable  invariant  properties  in  that  the 
constitutive  equations  and  the  equations  of  motion  are  properly  invariant 
under  arbitrary  (not  necessarily  infinitesimal)  superposed  rigid  body 
motions.  The  method  of  construction  is  effectively  such  that  they  consider 
the  rotation  at  one  material  point,  called  the  pivot  Y,  and  then  remove 
from  every  point  of  the  body  this  rotation  obtaining  a  new  configuration 
K  ,  and  then  show  that  in  the  configuration  <  all  constitutive  results 
are  properly  invariant.  In  the  second  part  of  their  work,  they  demonstrate 
tnat  to  the  order  of  approximation  of  the  kinematical  result,  it  does  not 
matter  which  particle  is  chosen  as  pivot. 

It  is  the  purpose  of  the  present  appendix  to  show  that  the  choice  of 

pivot  is  immaterial  also  in  the  presence  of  moderate  rotation  or  more 

specifically  for  (4.60),  (4.70)^,  (4.73)  and  (4.75).  Temporarily,  we 

.  *  ★ 

attach  the  subscripts  Y  and  Y'  to  quantities  F  ,H  ,  etc.,  when  the  pivots 
Y  and  Y '  are  identified  with  material  points  whose  position  vectors  are 
X  =  qX,  and  X=X  /  qX,  respectively,  in  the  reference  configuration  <Q. 

Using  (4.60)  and  (4.51),  the  rotation  tensor  at  X'  is  given  by 


Ry/.Rqc'.t)  =  R0(i  +  ^,  +  ;/)  . 

where 

=  4  (X',t)  =  0(eo)  as  -  0  , 

>F '  =  T  (X ',  t)  =  -  <F  /T  =  0  (e^)  ,  T  '  =  2$ '  . 

•w  ~  -V  0  — 

Recall  the  form  of  the  rotation  tensor  in  (4.60)  and  multiply  the  right- 
T 

hand  side  by  Ry /Ry /  =  I.  With  the  help  of  (Cl),  the  tensor  R  can  be 
written  as 


(Cl) 


(C2) 
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R  =  CRy/?y/^0Ci  +  t  +  P 


.  /  J.  Ml  t\ T  f 


5y/(I  +  !  + 1  >  (I  +  !  +  P 


=  Ry/Cl  +Q  +  l'  +  1-V'  ~l'l) 

=  r^{i+[*  +  $/-%(y'y  +  y  i')] 

+  [V  -  Y7  -kil'l  -  Y  Y7)]} 

3/2 

as  z.q  +  0,  where  terms  of  0(£q  )  as  eQ  ■*•  0  have  been  neglected.  Consider 


the  quantities 


*  =  *(X,t)  =  G  +  G'-lsCY'Y  +  Y  Y7)  =  <p  , 


Y  =  Y(X,t)  =  Y  -  Y7  -  %(Y  'Y  -  Y  Y')  =  -  Y 


^  =  .  wT 


and  observe  that  <t>  =  0 (e  1  and  Y*0(£s)  as  £  -*-0.  By  direct  calculation 
~  ~  o  *.  ~  o  o  ' 

3/2 

after  neglecting  terms  of  0(eQ  )  as  eo  +  0,  we  conclude  that 

Y  Y  *  2$  as  e  0  . 

~  ~  o 

In  addition,  both  4>  and  Y  vanish  at  X*X/.  Thus,  (C3)  has  the  form 


R  =  Ry  /Cl  +  $  +  Y)  » 

$  =  $T  =  0(e  )  as  e0  -*■  0  ,  Y  =  -  YT  =  0(ej  ,  ¥  Y  =  2¥  , 

Ry  /  =  R(X',t)  ,  f(X7,t)  =  Y(X',t)  ■  0  . 

Hence,  we  conclude  that  the  form  (4.60)  remains  unchanged  by  the  choice 
Y  ‘  as  the  pivot. 

From  (4.66)  and  with  the  help  of  (Cl)  and  (C4) ,  we  calculate 


- ir as.-. .»a AA. iti j  gjj.teWM1a,u,JateiA,,Jl  i,^„ ,  ..\,u  AA-JS, GU&Xu.  . -;:rg4ia"i 
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Zy'  =  V  ! 

=  (I  +  */  +  ¥/)TRTR  (I  +  <£>  +  f  +  E) 

..  ~  ~o~o  ~  ~  ~  ~ 

=  (I  +  $  +  $/  +  'F-'i,/-'i'/'F  +  E) 

=  I  +  $  +  ?  +E  ,  (C7) 

3/2 

where  terms  of  0(e  )  as  £  -*- C  have  been  neglected.  Thus, 

~  o  o 

H*z  =  F*/  -  I  =  T  +  Y+E  .  (C8) 

Using  CCS) ,  we  again  have 

Ey/  =  %(FyT  Fy/  -  I) 

=  k (By /  +  H* '  +  H H*  ,)  =  E  .  (C9) 

Hence,  (4.70) (4.73)  and  (4.75)  are  not  affected  by  the  choice  of 
pivot. 

In  the  above,  we  have  E  ^  =  0 (e  '“)  as  e  ->-0.  It  should  be  clear  that 

similar  procedures  can  be  used  for  the  case  in  which  E  =0(e  )  as  e  -*-0. 

^  ~,K  ~  o'  o 

Before  closing  this  appendix,  we  consider  the  consequence  of  the  change 

*  / 
cf  pivot  Y  on  the  stress  tensor  T  defined  by  (5.7).  When  a  particle  Y  , 

different  from  Y,  is  chosen  as  pivot,  (5.8)^  becomes 

T*,  =  FjJ,  T  Ry/  ,  (CIO) 

where  Ry  ,  =  R  (X=X/,t)  is  a  function  of  time  only.  Keeping  in  mind  that 
J([E]  -  0(eQ)  as  eo^0,  with  the  help  of  (Cl),  (C2)  ,  (3.7)  and  (S.8),  to  the 
order  of  approximation  considered  we  obtain 


Vhlr£ti&ifb  ibliiWiiiuiLu.  ^  j  ’<  /  •  J  f-lWJi.  W* 


Ceil) 


p  (I  +  «/  +  'F')RTR  JC  [E]rJr  (I  +$/  +  4'/) 
’o  ~  ~  -  ~o~o  ~  ‘■-.■‘-.o-o  ~  -  - 

P0JC  [E]  =■  T*  . 


3/2 

where  in  writing  (Cll)^  terms  of  0(eQ  )  as  eQ-*-Q  have  been  neglected. 

3/2  * 

Hence,  to  within  terms  of  0(eQ  ),  T  is  -maffected  by  a  change  of 
pivot. 
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